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Abstract 

Background: Accurate and precise estimates of glomerular filtration rate (GFR) are essential for clinical assessments, 
and many methods of estimation are available. We developed a radial basis function (RBF) network and assessed 
the performance of this method in the estimation of the GFRs of 207 patients with type-2 diabetes and CKD. 

Methods: Standard GFR (sGFR) was determined by 99m Tc-DTPA renal dynamic imaging and GFR was also estimated 
by the 6-variable MDRD eguation and the 4-variable MDRD equation. 

Results: Bland-Altman analysis indicated that estimates from the RBF network were more precise than those from 
the other two methods for some groups of patients. However, the median difference of RBF network estimates 
from sGFR was greater than those from the other two estimates, indicating greater bias. For patients with stage l/l I 
CKD, the median absolute difference of the RBF network estimate from sGFR was significantly lower, and the P 50 of 
the RBF network estimate (n = 56, 87.5%) was significantly higher than that of the MDRD-4 estimate (n = 49, 76.6%) 
(p < 0.0167), indicating that the RBF network estimate provided greater accuracy for these patients. 

Conclusions: In patients with type-2 diabetes mellitus, estimation of GFR by our RBF network provided better 
precision and accuracy for some groups of patients than the estimation by the traditional MDRD equations. 
However, the RBF network estimates of GFR tended to have greater bias and higher than those indicated by sGFR 
determined by 99m Tc-DTPA renal dynamic imaging. 
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Background 

Diabetic nephropathy is the leading cause of end stage 
renal disease, a condition characterized by abnormal 
glomerular filtration rate (GFR) and serum creatinine 
(SCr) [1]. The National Kidney Foundation (NKF) con- 
siders GFR as the best overall measure of kidney func- 
tion in healthy and diseased individuals [2]. However, 
measurement of GFR by use of radioisotopes is time- 
consuming and expensive, so this method is not used in 
routine clinical practice. Instead, numerous equations 
have been proposed to estimate the GFR without the 
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need for radioisotopes [2]. These equations consider SCr 
and several additional variables, such as age, gender, 
race, and body size [2]. The American Diabetes Associ- 
ation [3] also recommends estimation of GFR from 
serum creatinine (SCr) -based formulae, such as a 
Modification of Diet in Renal Disease (MDRD) study 
equation [4]. 

However, these equations may yield inaccurate esti- 
mates in some populations, such as elderly Chinese pa- 
tients with CKD [5]. Recent studies have criticized the 
equations currently used to estimate GFR in diabetic pa- 
tients [6-10]. In particular, the Chronic Kidney Disease 
Epidemiology Collaboration (CKD-EPI) equation [6], the 
Mayo Clinic Quadratic (MCQ) equation [7], and the 
four-variable MDRD equation [8] all underestimated 
GFR in patients with type-2 diabetes, and the Cockcroft- 
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Gault equation overestimated GFR in patients with type- 
2 diabetes [10]. These equations may be inaccurate 
because they do not account for ethnicity [11]. For 
example, in a group of Chinese patients with CKD, the 
MDRD equation 7 and the abbreviated MDRD equations 
underestimated GFR in patients with near-normal renal 
function and overestimated GFR in patients with ad- 
vanced renal failure [12]. These equations may also be 
inaccurate because they were developed by linear regres- 
sion methods [11,13,14]. Linear regression models do 
not account for the non-linear physiological processes 
that underlie GFR. Thus, it is important to develop bet- 
ter methodologies for estimation of GFR. 

Artificial neural networks (ANNs) have been success- 
fully used to model non-linear phenomena in the field of 
engineering forecasting. Modern ANNs provide effective 
nonlinear mapping of data, good fault tolerance, and 
good self-organization [15,16]. Previous research demon- 
strated that an ANN was more accurate than a logistic 
regression model in prediction of clinical outcome in pa- 
tients with systemic inflammatory response syndrome 
and hemodynamic shock [17]. Other research groups 
have used ANNs to estimate GFR, including a 
knowledge-based neural network (KBNN) model [18], 
an evolving connectionist systems (ECOS) model [19], 
and a tree-based model with 6 terminal nodes [15]. In 
all of these cases, the ANNs provided better estimates of 
GFR than the traditional equations. 

Radial basis function (RBF) networks are among the 
most widely used ANNs, but there have been limited 
clinical applications of these networks. Our previous 
study [20] described a simple RBF network for estima- 
tion of GFR (eGFR RBF ) in a group of 327 Chinese pa- 
tients with chronic kidney disease (CKD). The results 
indicated that the eGFR RB F had less bias and greater 
precision than the traditional MDRD equations. The 
accuracy (deviation less than 30% from the sGFR) of 
the eGFR RB F was significantly better than those from 
traditional eGFR equations, such as the Jelliffe-1973 
-equation and the Ruijin-equation [20]. 

In the present study, we tested the precision and ac- 
curacy of an RBF network model for estimation of GFR 
in an independent group of 207 Chinese patients who 
had type-2 diabetes and CKD and compared the results 
of the RBF network method with the results from two 
traditional MDRD formulae [13]. 

Methods 

Patients 

From January 2005 through December 2009, 207 con- 
secutive patients with type-2 diabetes from the Third Af- 
filiated Hospital of Sun Yat-sen University (Guangzhou, 
China) were enrolled. Patients younger than 18 years, 
taking cimetidine or trimethoprim, with acute kidney 



deterioration, clinical edema, skeletal muscle atrophy, 
pleural effusion or ascites, malnutrition, amputation, 
heart failure, or ketoacidosis were excluded. None of the 
patients were treated by dialysis during the study. CKD 
was staged according to the National Kidney Foundation 
(NKF) - Kidney Disease Outcomes Quality Initiative 
clinical practice guidelines [2] based on the GFR mea- 
sured by 99m Tc-DTPA dynamic imaging method. Pa- 
tients were placed into 3 groups based on CKD stage: (i) 
Stage I/II CKD (GFR > 60 mL/min/1.73 m 2 ); (ii) Stage 
III CKD (GFR = 30-59 mL/min/1.73 m 2 ); or (Hi) Stage 
IV/V CKD (GFR < 30 mL/min/1.73 m 2 ). The study 
protocol was approved by the institutional review board 
at the Third Affiliated Hospital of Sun Yat-sen University 
and all patients provided informed consent. All partici- 
pants provided written informed consent. 

Measurements 

GFR measured by the 99m Tc-DTPA renal dynamic im- 
aging method (modified Gate's method) was used as the 
standard GFR (sGFR) [21,22], and was calculated as 
described by Li et al. [23]. The gamma camera uptake 
method with 99m Tc-DTPA is a simple method for deter- 
mination of GFR, and is less time-consuming, and less 
expensive than other methods [24]. Moreover, this 
method has been recommended as the reference ap- 
proach for determination of GFR by the Nephrology 
Committee of Society of Nuclear Medicine [25], and is 
widely used as a standard method for evaluation of kid- 
ney function and estimation of GFR in China. m Tc- 
DTPA renal dynamic imaging was measured by a Mil- 
lennium TMMPR SPECT using the General Electric 
Medical System, as described previously [22]. Serum al- 
bumin (Alb) and blood urea nitrogen (BUN) were 
assayed on a Hitachi 7180 autoanalyzer (Hitachi, Tokyo, 
Japan; reagents from Roche Diagnostics, Mannheim, 
Germany). SCr was measured by an enzymatic method 
on the same instrument according to the manufacturer's 
instructions. SRM967 (standard reference material re- 
leased by NIST for serum creatinine calibration) was 
used for calibration. Patient sex, age, height, and weight 
were recorded at the same time. 

RBF network 

An ANN is a computational method composed of 
interconnected artificial neurons (mathematical func- 
tions) that processes information and that typically con- 
sists of an input layer, one or more hidden layers, and an 
output layer. ANNs are used in diverse scientific and en- 
gineering fields to model the complex relationships of 
inputs and outputs. An RBF network is a feed-forward 
network with one hidden layer, in which activation of 
the hidden layer is a nonlinear radial basis function 



Liu et al. BMC Nephrology 201 3, 14:1 81 
http://www.biomedcentral.com/1471-2369/14/181 



Page 3 of 9 



(a function whose value only depends on the distance to 
the origin). 

In this study, the input layer consisted of measured 
serum creatinine (SCr) and the output layer consisted of 
sGFR. Our previous work [20] indicated that when SCr 
was measured by the enzymatic method, a simple RBF 
network model successfully estimated GFR in a popula- 
tion of 327 Chinese patients with CKD, based on ana- 
lysis of all patients and on analysis of subgroups of 
patients with different stages of CKD [20]. In this previ- 
ous study, the RBF network was a feed-forward ANN 
with an input layer of one unit (SCr), one hidden layer, 
and an output layer of one unit (sGFR), which was mea- 
sured in all 327 CKD patients. This RBF network was 
constructed by use of the newrbe function in 
Math Works. In the present study, we tested this RBF 
network in an independent group of 207 patients who 
had type-2 diabetes (external validation data set) to ver- 
ify the original results. 

MDRD equations 

GFR was also estimated by the traditional MDRD equa- 
tions [13]. In particular, we used the re-expressed 4- 
variable MDRD equation (R-MDRD4): 

eGFR 4 = 175 x Scr" 1154 x Age 0203 

x [0.742 if female] x [1.212 if black] 

and the re-expressed 6-variable MDRD equation 
(R-MDRD6): 

eGFR 6 = 161.5 x Scr 0 " 9 x Age" 0176 x BUN" 017 
x Alb +0 ' 318 x [0.762 if female] 
x [1.18 if black] 

Statistical analysis 

All demographic and clinical data are summarized as 
means ± standard deviations (SDs), as medians and 
inter-quartile ranges (IQRs: Ql, Q3) for continuous vari- 
ables that had non-normal distributions, and as N and 
percent for categorical data (CKD stage). Data were 
compared using a one-way ANOVA with Bonferroni's 
post hoc correction, the Kruskall Wallis test, the Mann- 
Whitney U test for pair-wise comparisons of data that 
had non-normal distributions, or Pearson's Chi-square 
test (gender). 

The overall differences between eGFR and sGFR are 
summarized as medians and IQRs due to the non- 
normal distributions. Differences among patients with 
different CKD stages were compared with the Kruskall 
Wallis test with a post hoc method or with the Mann- 
Whitney U test for pair-wise comparisons. Within-group 
comparisons of measurements were performed using the 
Wilcoxon signed ranks test for a given CKD stage. 



The accuracy of eGFR are summarized as N and per- 
cent of patients with eGFR differing less than 15% (P15), 
30% (P 30 ), and 50% (P 50 ) from sGFR. Accuracy of the es- 
timates was compared for patients with different CKD 
stages using the Pearson Chi-square test. The accuracies 
of eGFR values were compared using the McNemar test 
within the same CKD stage. Bland- Altman plots (eGFR 4 
vs. sGFR, eGFR 6 vs. sGFR, and eGFR RB F vs. sGFR) were 
graphed with Medcalc for Windows (ver. 9.3.9.0, 
Mariekerke, Belgium). The 3 different methods of esti- 
mating GFR were also used to classify patients by CKD 
stage. A Wilcoxon sign-rank test was used to compare 
the differences of CKD stages from sGFR and each of 
these estimates. 

All statistical assessments were two-tailed and a p-vahie 
less than 0.05 was considered significant. The significance 
level was adjusted by Bonferroni's method to 0.0167 (0.05/ 
3) and 0.01 (0.05/4) for post hoc pair-wise comparisons of 
CKD stages and eGFR, respectively. All statistical analysis 
was performed using SPSS (version 11.0 SPSS, Chicago IL, 
USA). 

Results 

A total of 207 patients (119 males, 88 females) with a 
mean age of 61.43 years (SD = 12.03) were enrolled. 
Table 1 shows the baseline demographic and clinical 
characteristics of the 207 patients and of sub-groups 
with different stages of CKD. There were 64 patients 
(30.9%) with stage I/II disease, 81 patients (39.1%) with 
stage III disease, and 62 patients (30%) with stage IV/V 
disease. As expected, patients with more severe disease 
had lower serum Alb and sGFR, and higher SCr and 
BUN (Table 1). 

Table 2 shows the sGFR values and the three different 
estimates of GFR (eGFR 4 , eGFR 6 , and eGFR RBF ) for all 
207 patients and for sub-groups with different stages of 
CKD. As expected, GFR values were lower in patients 
with more advanced disease. The eGFR RB F of all 207 
CKD patients and of patients with different stages of 
CKD were significantly higher than the sGFR values (p < 
0.01 for all). In addition, for patients with stage III or 
stage IV/V CKD, the eGFR RBF values were significantly 
higher than those from eGFR 4 and eGFR 6 (p < 0.01 for 
both). 

The validity of an estimate is a function of its accuracy 
and precision, and a valid estimate should be close to 
the true value (low bias) and be reproducible [26]. 
Table 3 presents the overall performances of the 3 differ- 
ent estimates of GFR in all 207 patients and in patients 
with different stages of CKD. Each row in the table 
shows the median difference of eGFR and sGFR and the 
inter-quartile range (IQR), the median absolute differ- 
ence of eGFR and sGFR and the IQR, and the percent of 
GFR estimates within 15% (P 15 ), 30% (P 30 ), and 50% 
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Table 1 Demographic and clinical characteristics of patients with type-2 diabetes and different stages of CKD (n 


= 207) 


Variable 


Total 1 (n = 207) 


Stage l/ll 1 (n = 64) 


Stage lll 1 (n = 81) 


Stage IV/V 1 (n = 62) 


P-value 


Age (years) 


61.43 ± 12.03 


56.64 ± 1 1 .65 


62.78 ± 11. 53 a 


64.61 ±11.72 a 


<0.001* 


Sex 










0.177 


Males 


119 (57.5) 


34 (53.1) 


53 (65.4) 


32 (51.6) 




Females 


88 (42.5) 


30 (46.9) 


28 (34.6) 


30 (48.4) 




BMI (kg/m 2 ) 


23.42 (2145 , 25.83) 


23.54 (21.32 , 25.77) 


22.86 (21.16 , 25.52) 


24.05 (22.04 , 25.92) 


0.262 


BSA <m 2 ) 


1 .66 (1 .54 , 1 .80) 


1.66 (1.52 , 1.77) 


1.64(1.54, 1.81) 


1.69 (1.54 , 1.77) 


0.826 


Alb (g/dL) 


3.75 (3.19,4.20) 


4.00 (3.51 , 4.41) 


3.83 (3.30 , 4.18) 


339 (2.99 , 3.84) ab 


<0.001* 


SCr (mg/dL) 


1.52 (0.92 , 3.72) 


0.77 (0.57 , 1 .09) 


1.58 (1.09 , 2.39) a 


5.19 (2.99 , 7.01) ab 


<0.001* 


BUN (mg/dL) 


27.17 (17.76,48.54) 


16.70 (13.74 , 21.25) 


27.90 (18.74 , 44.82) a 


54.68 (41.29 , 70.60) ab 


<0.001* 


DTPA-GFR (mL/min/1.73 m 2 ) 


43.45 (2535 , 64.36) 


78.39 (66.86 , 87.89) 


43.42 (35.81 , 54.82) a 


19.70 (14.58 , 23.1 9) ab 


<0.001* 


GFR §60 mL/min/1.73 m 2 


64 (30.9) 


64 (100) 


0(0) 


0(0) 




GFR <60 mL/min/1.73 m 2 


143 (69.1) 


0 (0) 


81 (100) 


62 (100) 





Abbreviations: BMI body mass index, BSA body surface area, Alb albumin, SCr serum creatinine, BUN blood urea nitrogen, DTPA-GFR GFR measured by m Tc-DTPA 
renal dynamic imaging. 

1 Continuous data are given as means ± SDs for age, as medians and IQRs for continuous variables with non-normal distributions, and as n (%) for categorical 
variables. Differences among patients with different stages of CKD were compared with a one-way ANOVA and a Bonferroni post-hoc comparison for age, a 
Kruskall Wallis test with Mann-Whitney U test for pair-wise comparison of other not normally distributed continuous data; and Pearson's Chi-square test for sex. 
*p < 0.05, indicates significant difference among CKD stages. 

a b p< 0.0167 (0.05/3), indicates significant difference of CKD stage l/ll and III, respectively. 



(P 50 ) of sGFR [26]. The results indicate that for all 207 
patients, the IQR of the eGFR RBF was smaller than those 
from eGFR 4 and eGFR 6 , indicating better precision for 
the RBF network; however, the median difference of the 
eGFR RBF from sGFR was significantly higher than those 
from eGFR 4 and eGFR 6 (p < 0.0167 for both). The same 
trends occurred in patients with stage III CKD and stage 
IV/V CKD. The eGFR RBF had a lower median absolute 
difference from sGFR and a higher percent of estimates 
within 15% (P 15 ) and 30% (P 30 ) of sGFR, suggesting that 
eGFR RBF had better accuracy. However, the differences 
between eGFR RBF and eGFR 4 and eGFR 6 were not statis- 
tically different. The results also indicate that for the pa- 
tients with stage I/II CKD, the median absolute 
difference of eGFR RBF and sGFR (19.73) was significantly 
less than that of eGFR 4 and sGFR (27.28) (p < 0.0167), 
although the P 50 for eGFR RBF (n = 56, 87.5%) was signifi- 
cantly higher than that for eGFR 4 (n = 49, 76.6%) (p < 
0.0167). 

Figure 1 shows Bland- Altman plots for comparisons of 
sGFR with GFR estimated by eGFR 4 (A), eGFR 6 (B), and 
eGFR RBF (C). The precision is indicated by the distance 
between the dashed lines (95% limits of agreement) [27]. 
The results indicate that the distance between lines of 
95% limits of agreement was 71.4 for eGFR RBF) signifi- 
cantly lower than that for eGFR 4 (91.5) and eGFR 6 
(85.2). Thus, this analysis also indicates that the eGFR RBF 
had better precision than the two other estimates from 
the MDRD equations. However, Figure 1 also shows that 
the mean differences of sGFR and eGFR 4 was 3.1, sGFR 
and eGFR 6 was 1.0, and sGFR and eGFR RBF was 11.1. 
This indicates that eGFR RBF had greater bias than the 



estimates from the MDRD equations and over-estimate 
the GFR. 

Discussion 

We compared the performance of an RBF neural net- 
work in the estimation of GFR with the performance of 
two traditional GFR estimates based on the MDRD 
equations (MDRD-4 and MDRD-6) in patients with 
type-2 diabetes and different stages of CKD. Our results 
indicate that the RBF network provided more precise es- 
timates of GFR than the MDRD equations, and also pro- 
vided significantly more accurate estimates of GFR for 
patients with stage I/II CKD. However, the RBF network 
also had higher bias than the traditional MDRD equa- 
tions. In particular, the eGFR RBF tended to over-estimate 
GFR more than eGFR 4 and eGFR 6 , especially for patients 
with CKD stage IV/V (Table 2). 

In the field of medical data processing, the theoretical 
basis for the use of statistical regression methods is the 
"law of large numbers". That is, the difference of the 
average of many measurements from the true value 
should be smaller as more measurements are recorded. 
However, application of a model derived from one data 
set to another data set may yield poorer accuracy. More- 
over, regression methods can only be used for a limited 
number of models, and interactions among variables 
places limits on their use. ANNs have no a priori re- 
quirement for data distribution, and can handle multi- 
collinear input variables, neither of which can be man- 
aged by regression methods. 

These advantages of ANNs have led to their use in 
several previous studies for estimation of GFR. Song 



Table 2 Glomerular filtration rates (sGFR, eGFR 4 , eGFR 6 , and eGFR RBF ) of patients with different stages of CKD 



Measurement method (mL/min/1.73 m 2 ) 


Total 1 (n = 207) 


Stage l/ll 1 (n 


= 64) 


Stage lll'fn 


= 81) 


Stage IV/V 1 (n = 62) 


P-value 


sGFR 


43.45 (25.35 , 64.36) 


78.39 (66.86 , 


87.89) 


43.42 (35.81 


, 54.82) a 


19.70 (14.58 , 23.1 9) ab 


<0.001* 


eGFR 4 


42.43 (15.68 , 75.46) 


89.12 (67.73 , 


1 24.20) + 


42.43 (26.50 


, 59.31 ) a 


9.60 (6.71 , 17.11) ab * 


<0.001* 


eGFR 6 


39.98 (15.28 , 72.28)* 


90.88 (64.52 , 


116.84)* 


40.21 (24.72 


, 60.54) a * 


9.92 (6.54 , 1 7.02) ab * 


<0.001* 


eGFR RBF 


52.25 (34.21 , 81.64) +M 


94.23 (70.34 , 


115.27)* 


50.42 (39.61 


, 70.57) a ** 5 


24.01 (18.70 , 37.28) ab ** § 


<0.001* 



Abbreviations here and in Tables 3: sGFR standard GFR, eGFR 4 GFR estimated by re-expressed MDRD-4, eGFR 6 GFR estimated by re-expressed MDRD-6, eGFR RBF GFR estimated by RBF network. 

1 Data are given as medians and IQRs for continuous variables because of non-normal distributions; differences among CKD stages were compared using the Kruskall Wallis test with a post hoc method, or the Mann- 
Whitney U test for pair-wise comparisons; within-group comparisons of the four measurements were performed using Wilcoxon Signed Ranks Test for a given CKD stage. 
*p < 0.05, indicates significant difference among CKD stages. 

a,b p < 0.01 67 (0.05/3), indicates significant difference in a comparison of CKD stage l/ll and III, respectively. 
+, * ,§ , p < 0.01 (0.05/4), indicates significant difference in a comparison with sGFR, eGFR 4 , and eGFR 6 , respectively. 



Table 3 Overall performance of different methods used to estimate glomerular filtration rate in patients with different stages of CKD 





Median difference of eGFR and sGFR^ 


Median absolute difference of eGFRand sGFR 2 


15% (P-.il 


30% (Pan) 
/u \r 30/ 


50% (Pen) 


Total (n = 207) 














eGFR 4 and sGFR 1 


-2.83 (-11.23 , 12.13) 


32.54 (15.74 , 


58.97) 


49 (23.7%) 


92 (44.4%) 


140 (67.6%) 


eGFR 6 and sGFR' 


-4.06 (-12.59 , 1 1.21)t 


32.78 (16.43 , 


57.38) 


49 (23.7%) 


91 (44.0%) 


139 (67.1%) 


eGFR RBF and sGFR 


9.76 (-0.19 , 18.02)++ 


26.24 (13.29 


55.82) 


62 (30.0%) 


112 (54.1%) 


149 (72.0%) 


Stage l/ll (n = 64) 














eGFR 4 and sGFR 1 


9.23 (-9.20 , 37.35) 


26.41 (12.20 


49.65) 


1 9 (29.7%) 


35 (54.7%) 


49 (76.6%) 


eGFR 6 and sGFR 1 


4.13 (-14.66 , 32.62)t 


27.28 (12.73 , 


44.42) 


1 6 (25.0%) 


33 (51.6%) 


49 (76.6%) 


eGFR RBF and sGFR 


1 2.03 (-4.98 , 27.83) 


19.73 (12.98 , 


35.56)+ 


21 (32.8%) 


43 (67.2%) 


56 (87.5%)+ 


Stage III (n = 81) 














eGFR 4 and sGFR 1 


-1.04 (-12.85 , 9.92) 


27.01 (13.76 


. 58.47) 


22 (27.2%) 


43 (53.1%) 


59 (72.8%) 


eGFR 6 and sGFR 1 


-4.80 (-14.07 , 8.97)+ 


29.82 (13.65 


57.57) 


25 (30.9%) 


41 (50.6%) 


58 (71.6%) 


eGFR RBF and sGFR 


7.96 (-1.22 , 21.47)++ 


25.40 (11.23 


48.50) 


26 (32.1%) 


47 (58.0%) 


61 (75.3%) 


Stage IV/V (n = 62) 














eGFR 4 and sGFR 1 


-6.48 (-11.06,-0.21) 


46.90 (30.81 , 


64.90) 


8 (12.9%) 


14 (22.6%) 


32 (51.6%) 


eGFR 6 and sGFR 1 


-6.99 (-11.23,-1.63) 


46.49 (26.56 , 


65.17) 


8 (12.9%) 


1 7 (27.4%) 


32 (51.6%) 


eGFR RB F and sGFR 


8.89 (1 .95 , 1 5.48)++ 


49.20 (14.98 


89.37) 


1 5 (24.2%) 


22 (35.5%) 


32 (51.6%) 



1 Abbreviations as in Table 2. 

2 Data are given as median difference (IQRs) and median absolute difference (IQRs) between eGFR and three types of eGFR. Within-group comparisons between the measurements were performed using the Wilcoxon 
Signed Ranks Test for a given CKD stage. 

3 The accuracy of eGFR was expressed as n (%) of measurements within 15% (Pis), 30% (P30), and 50% (P 50 ) of sGFR. Data within-measurements were compared using the McNemar test for a given CKD stage. 
*p < 0.05, indicates significant difference among CKD stages. 

+ '*, P< 0.01 67 (0.05/3), indicates significant difference with eGFR 4 and sGFR and eGFR 6 and sGFR, respectively. 
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Figure 1 Bland-Altman plots of the estimation of GFR (mL/min/1.73 m 2 ) relative to sGFR by the (A) re-expressed 4-variable MDRD 
equation (eGFR 4 ), (B) re-expressed 6-variable MDRD equation (eGFR 6 ), and (C) RBF network (eGFR RB F)- Solid lines represent mean 
differences and dashed lines represent 95% limits of agreement of the mean difference. 



et al. [18] used a knowledge-based neural network model 
(KBNN) for evaluation of renal function based on 441 
GFR data vectors from 141 patients. Their proposed 
GFR prediction model had at least 10% better accuracy 
than any of the individual regression formulae or a 
standard neural network model. Marshall et al. [19] used 
evolving connectionist systems (ECOS), in which com- 
puting structures are trained to generate output from a 
given set of input variables. They concluded that their 
ECOS model provided better prediction of GFR in rou- 
tine clinical practice. No ANNs have been used to esti- 
mate GFR of patients with type-2 diabetes mellitus. 

It is noteworthy that the RBF network used in this 
study was only based on SCr measurements, in contrast 
to the re-expressed MDRD equations, which require 
measurement of 4 or 6 variables. The NKF does not rec- 
ommend use of SCr alone for assessment of kidney 
function [2]. However, previous research indicated that 
use of fewer variables can yield acceptable estimates of 
GFR. For example, Bevc et al. [28] reported that a 
cystatin C-based estimate, which only requires measure- 
ment of serum cystatin C, is a reliable marker of GFR in 
elderly patients and is comparable to the creatinine- 
based formulae, including the CKD-EPI formulae. Our 
results suggest that an RBF model based on a single 
measurement (SCr) can provide precise and accurate es- 
timates of GFR. 



There are several limitations to our study. First, SCr 
was measured by the enzymatic method. Peake et al. 
[29] indicated that the enzymatic creatinine assay, al- 
though theoretically more specific, can have interference 
problems. However, this method produces results for pa- 
tient samples that agree closely with the results from 
isotope dilution mass spectrometry (ID-MS). This moti- 
vated our use of the re-expressed MDRD equations 
(MDRD-4 and MDRD-6) instead of the original equation 
[4], because the original MDRD equation was developed 
for use with ID-MS traceable serum creatinine [30]. Sec- 
ond, ANN models can be difficult to display as equa- 
tions and cannot be readily used without special 
software, so physicians may be reluctant to accept the 
results of ANN models. There is need for a platform that 
can display ANN models and that allows other re- 
searchers to readily perform external validation. Third, a 
previous study-indicated that GFR estimated by 99mTc- 
DTPA dynamic renal imaging might not better than the 
modified abbreviated MDRD equation [31], and the 
renal dynamic imaging method was less accurate than 
CKD-EPI equation as well [32]. However, the same study 
found that the two methods performed similar capability 
in determining GFR among higher-GFR patients [32], 
and 99m Tc-DTPA dynamic renal dynamic imaging yields 
accurate results that are nearly the same as those from 
measurements of inulin clearance [33]. Rehling et al. 
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showed that a regression line between the values mea- 
sured by these different methods did not differ from the 
line of identity [33]. Ma et al. (2007) suggested that, 
using proper reference GFR, more adequate background 
subtraction, and soft-tissue attenuation correction may 
improve the accuracy of 99m Tc-DTPA dynamic renal im- 
aging [31]. Finally, our RBF network predicted a higher 
GFR than that from 99m Tc- DTPA renal dynamic im- 
aging. This might be due to differences of participants in 
the training group (CKD patients with and without dia- 
betes [20]) and the study group (diabetes patients with 
and without normal kidney function). Use of more simi- 
lar training and study groups would provide better exter- 
nal validation and may provide improved results. 

A recent survey in China [34] showed that the preva- 
lence of diabetes was 9.7%, corresponding to 92.4 million 
people. Although some of the established methods used 
to estimate GFR are suitable for Chinese patients with 
CKD [22], it is important to have more accurate and 
precise estimations of GFR. In some measures of accur- 
acy and precision, our RBF neural network performed 
significantly better than the re-expressed MDRD equa- 
tions in the estimation of GFR. In particular, the IQRs 
(Table 3) and 95% limits of agreement (Figure 1) for the 
eGFR RB F were smaller than those from eGFR 4 and 
eGFR 6 , indicating better precision for the RBF network. 
However, our data indicated that eGFR estimated by the 
RBF neural network tended to be higher than the sGFR, 
and this would result in under-estimation of CKD stage. 
We suggest that use of an RBF network model with 
more variables and testing of the model with additional 
data sets may ultimately provide more accurate and pre- 
cise estimates of GFR. 

Conclusions 

In patients with type-2 diabetes, GFR estimated by our 
RBF network provided better precision and accuracy for 
some groups of patients than GFR estimated by the trad- 
itional MDRD equations. However, the RBF network es- 
timates of GFR tended to have greater bias and higher 
than those indicated by sGFR determined by 99m Tc- 
DTPA renal dynamic imaging. 
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